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In this paper, we provide a theoretical analysis of strongly interacting quantum systems confined 
by a time-dependent external potential in one spatial dimension. We show that such systems can be 
used to simulate spin chains described by Heisenberg Hamiltonians in which the exchange coupling 
constants can be manipulated by time-dependent driving of the shape of the external confinement. 

As illustrative examples, we consider a harmonic trapping potential with a variable frequency and 
an inhnite square well potential with a time-dependent barrier in the middle. 

PACS numbers: 67.85.-d 75.10.Pq 


Introduction. ~ Strongly interacting quantum systems 
are an intricate and exciting part of theoretical physics. 
Their intricacy is due to the strong many-body correla¬ 
tions that may lead to unexpected new phenomena not 
present in the weakly-interacting case. 

For systems with strong interparticle coupliM, one 
spatial dimension (ID) plays a very special role [l|. One 
reason for this is the unusual duality, often called the 
Fermi-Bose mapping, between ID impenetrable bosons 
and ideal fermions, which was rigorously shown in 1960 
by Girardeau [H. The most exciting aspect of this du¬ 
ality is the possibility to study it in modern experimen¬ 
tal setups with two different atomic species As a 

future perspective, the Fermi-Bose mapping suggests 
to engineer a chain of spins with adjustable nearest- 
neighbor couplings using a strongly repulsive multicom¬ 
ponent system in a trap [iMl- Such spin chains pos¬ 
sess a very high degree of tunability thus opening the 
possibility of realizing and studying phenomena such as 
ID SUJN) quantum magnets and perfect state trans¬ 
fer @,[11 . 


While the Fermi-Bose mapping was first established 
for a stationary system, the generalization to the case 
of a time-dependent trapping potential is straighforward 
for a system of impenetrable bosons [3, [3l- To the 
best of our knowledge, however, such a generalization for 
multicomponent systems with large but finite interaction 
was not previously discussed in the literature [3|- Due 
to the interplay of two different time-dependent effects, 
this generalization is far from obvious. First, there is the 
motion of particles due to the time-dependent trapping 
potential, and second, there is the particle exchange. As 
we will show the timescales for these effects are effectively 
decoupled from one another and the dynamics of particle 
exchange is determined by the trapping potential. 

In this paper, we consider a system with two kinds 
of spinless fermions with strong interspecies repulsion. 
We first show that the behaviour of such a system can 
be described by the Heisenberg Hamiltonian with time- 
dependent exchange coupling coefficients. These coeffi¬ 


cients can be altered by manipulating the shape of the 
trapping potential as a whole. This contrasts our idea 
with an idea of realizing a time-dependent Heisenberg 
Hamiltonian on a lattice by addressing every site inde¬ 
pendently. As we discuss below, our approach has very 
different strengths and limitations and thus ideally com¬ 
plements the standard lattice approach. In particular, it 
allows to address any trapping potential and is not lim¬ 
ited to the lattice approximation. For a four-atom system 
this is sketched in Fig.[T] Part a) shows the initial config¬ 
uration with two fermions in one well and two fermions of 
a different kind in the other well. We also sketch a possi¬ 
ble evolution of this configuration in a time-dependent 
potential where the final state corresponds to the ex¬ 
change of the two pairs. This evolution can be described 
by mapping the system to a spin-chain model described 
by a Heisenberg Hamiltonian where the coupling coeffi¬ 
cients depend on time, see Fig. □b). The mapping opens 
a way to engineer and simulate driven Heisenberg Hamil¬ 
tonians with time-dependent coefficients where the time 
dependence gives an extra knob to tune the dynamics in 
the system 17. 3. 

To illustrate our findings, we apply this mapping to a 
time-dependent harmonic potential and an infinite square 
well potential with a time-dependent barrier in the mid¬ 
dle. For the former case, we show that the coupling co¬ 
efficients in the spin chain are simply multiplied with a 
position independent scale factor. In the latter case, one 
can tune the middle coupling coefficient almost indepen¬ 
dently from the others. This allows one to achieve a 
controlled exchange of pairs, see Fig. [T]a). 

Formulation. - For the sake of the argument let us 
start with a ID system of N spinless fermions of one 
kind (spin up) and one fermion of another kind (spin 
down). We assume that every particle has mass m and is 
confined by the same time-dependent trapping potential 
eV{x/L,et/h), where e = /(jnL'^) and L is some natu¬ 

ral time-independent unit of length. For convenience, we 
assume that m = h = L = 1 from now on. 

The dynamics of a system with a spin-up fermion 
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FIG. 1. (Color Online). Mapping of a strongly interacting 
one-dimensional system in a time-dependent potential , a), 
onto a spin chain with time dependent coefficients, b). 


placed at x and spin-down fermions at yi,...,yN is de¬ 
scribed by the wave function t/i,..., j/at, t), which sat¬ 
isfies the Schrodinger equation, 


dt 


N 

H = '^h{yi,t) + h{x,t), ( 1 ) 


where h{x,t) = + V{x,t) is the one-body Hamil¬ 

tonian. The zero-range interaction enters through ID 
Bethe-Peierls boundary conditions at the points where 
the particles meet (see e.g. Ref. 0 ): 




/d'i> 

VSa; dyj - 


2g'i>{x = yi), 


( 2 ) 

where g is the interaction strength and the notation 
X = yf^ means that the derivative is taken at the point 
X = yi zk e, with e > 0 and the limit e ^ 0 is taken 
afterwards. Below, we consider the dynamics of the sys¬ 
tem in the following scenario: the interaction is adiabati- 
cally tuned in a constant trapping potential V {x, 0) from 
zero to some value g^ which is very large. This proce¬ 
dure initializes the state j/i,..., ^at). It is assumed 

that at later times the shape of the trapping potential 
depends on time and we look for the wave function satis¬ 
fying Eq. O with the initial condition l['(t = 0 ) = 


Initial state. - Let us start by discussing the initial 
wave function If 1/g^ = 0 then Eq. ([2]) dictates 

that the particles cannot exchange their relative positions 
and should be described separately on each ordering 
of particles, e.g. x < yi < y 2 < ■■■ < yN, on which the 


solution is obtained from the Fermi-Bose mapping 


i: 


AT+l 

^( 0 ) ^ '^af\j{x,yi,...,yN)<k^°Hx,yi,...,yN), (3) 


i=i 


where function Xj is non-zero only if it contains j — 1 
arguments yi that are smaller than x, and is one 
of the eigenstates of the Hamiltonian for -|- 1 spinless 
fermions (for the illustrative examples below we use the 
ground state). First note that the states from Eq. ([3]) 
are N + 1-fold degenerate even if the eigenspectrum of 
spinless fermions is non-degenerate. Thus, to find '1'°, 
we should find the adiabatic eigenstates in g, that are 
characterized by This can be done perturbati yely 
by minimizing the energy in the limit 5 / —>■ oo [i 0 ,[ 2 l[. 
For large but finite interaction strengths, the wave func¬ 
tion preserves the form given by Eq. © but acquires an 
additional contribution proportional to 1/g/. Further¬ 
more, the minimization of energy leads to the mapping 
of a system onto a spin chain. To establish such a map¬ 
ping in the time-dependent case, where the energy is not 
a good quantum number, a new approach is necessary 
and this is what we provide in this paper. 


Time dynamics. - At t > 0, the external potential 
depends on time and the time evolution is described by 
'l'(t). Let us first consider the system with infinite in¬ 
teraction, i.e. 1 / 5 / = 0. In this case, the wave func¬ 
tion at each ordering should still be described with the 
wave function of spinless fermions 0 , <k{x,yi,...,yN,t). 
Moreover, the probability of each ordering cannot be 
changed since the particles do not exchange their posi¬ 
tion. So d' (t) in this limit has the same form as in Eq. (|3l) 
with $(t) instead of 

Let us now assume that the interaction strength is 
large but finite. Apparently this means that the wave 
function at each ordering cannot be described exactly 
with and we should to look for a solution in the 

form = cj) + —f where the function cj) reads 


AT-1-1 


cj){x,yi,...,yN,t) = ^ aj{t)xj{x,yi,...,yN)<k. (4) 


1=1 


Without any loss of generality, we assume that 
/ /*Xi‘J’d 2 :dyi...dyAf = 0 , i.e. that the functions / and 
4) are orthogonal on each ordering of the coordinates. 
Having in mind these conditions, we insert Eq. o in 
the Schrodinger equation. To proceed further, we insert 
the ansatz wave function into the Schrodinger equation 
and project it on each ordering. Next using that / is 
orthogonal to <f> at every moment of time together with 
the boundary conditions (I 2 |), we eliminate the function 
/ (See Supplemental Material 22|). This procedure al¬ 
lows us to obtain a system of equations for the set of 
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coefficients, aj(t), 


i — aj{Jj-i + Jj) aj-iJj-i Qj^iJj +0 ^^ 2 ^ ’ 

^ (5) 

where, assuming that Jq = Jn+i = 0 , the parameters Ji 
are defined as follows 


Mt) = - 


if. 


yi<y2<-<yN 


dyi-.AyN 


d<s>{t) 

dx 


2 

x^Vj 


9f L 


x<yi<V2<---<yN 


dxdyi...d?/Ar$2(i) ' 


( 6 ) 


After writing Eq. ([5]) in matrix form, it becomes appar¬ 
ent that up to the order 1 /this equation also describes 
the dynamics of a spin chain with the Heisenberg Hamil¬ 
tonian 


= (7) 

t=i 

and the corresponding wave function is 
v-i-i 

-tiv+i), ( 8 ) 


where we denote the identity operator on every site 
with I, cr-J = {al,(T^,ai) are the Pauli matrices acting 
on a spin at site j, and Jj are site- and time- dependent 
interaction coefficients. Equations (0 and ([H]) generalize 
the time-independent mapping [ 6 l-l 8 | onto a spin-chain 
Hamiltonian to the time-dependent case. The derivation 
above implies that the time scale for the particle mo¬ 
tion in leading order (in 1/gf) is determined by the trap 
alone, whereas the time scale of the spin exchange is pro¬ 
portional to 1/qf. It is related to the famous spin-charge 
separation |23l| in ID, although here we derived it for a 
strongly interacting mesoscopic system from first princi¬ 
ples in the presence of an external potential that depends 
on time. Note that there are higher order contributions 
to both, the particle motion and spin exchange. How¬ 
ever, these corrections are negligible in the case of strong 
interactions, 1/gf 0 , and therefore we do not need to 
consider them here. 

Applying the presented approach it is easy to show that 
the Hamiltonian Q , can be used for any number of spin- 
down fermions similar to the time-independent case (see 
Ref. @ for a derivation). This is due to the fact that the 
main process in the system is the spin exchange of neigh¬ 
boring particles which is correctly described in the Hamil¬ 
tonian 0. The same logic also applies to multicom¬ 
ponent system or systems made of strongly-interacting 
bosons. 

Discussion. - We first assume that the coupling coef¬ 
ficients, Ji, are independent of time. Then linear system 
of equations ([5]) has the fundemental set of solutions: 
a jit) = where e is the relevant eigenvalue of 


the Hamiltonian ([7]). Let us now consider what happens 
if the external trapping potential depends on time. To 
find the coefficients Ji in this case, we first need to solve a 
time-dependent one-body problem and construct a Slater 
determinant wave function $ out of the established solu¬ 
tions. 

As our first application, we consider a system trapped 
by a harmonic oscillator potential, H(a:,D = u;^(t)x^/ 2 , 
for which one-body solutions are known [241427l | , yielding 
^ix,yi,...,yN,t) from as 




$ = 




\/Ht) 


AT-1-1 







(9) 


where E is the initial energy and A(t) is the time- 
dependent scale parameter. Its time derivative A is 
determined from the equation: A^A = I — uj'^ (t) (t). 

Since our choice of units sets u;^(t < 0) = I, the ini¬ 
tial conditions for this equation read X{t < 0) = 1 and 
A(t < 0) = 0. Obviously, such a wave function $ pro¬ 
duces Ji(t) I/A^(t), so that all coupling constants de¬ 
pend on time in the same way. The corresponding system 
of equations ([5]) has the following fundamental set of so¬ 
lutions: aj{t) = cL'^^e Wo. Thus we see that the 
scale invariance given by the harmonic trap is preserved 
up to terms suppressed by 1/gf in the form: 


T = 






Therefore, the overall spin dynamics in the system is not 
affected by a change of the external potential up to cor¬ 
rections suppressed hy 1/gf. Of course, the harmonic os¬ 
cillator is a truly special case due to the scale invariance 
and any trapping potential that is not scale invariant will 
have more pronounced effects on the system. A detailed 
discussion of the breaking of scale invariance in the os¬ 
cillator for two particles by higher order corrections can 
be found in Ref. 281. 

It is interesting to note that the spin dynamics of the 
system in a harmonic trap can be altered by a time- 
independent weak magnetic field Q. With a magnetic 
field the Hamiltonian is H = H + 


9f 


where, for simplicity, we again consider a system with 
only one spin-down fermion [ 2 ^ . The correponding spin 
chain Hamiltonian is written as 


1 ^ 


AT+l 

E 

i=i 




( 11 ) 


where 

3 it) - ^ /d2/i-dyivda:Xj-(x,?/i,...,yAf)|$p&(a:) 

^ 9f J dyi...dyNdxXjix,yi,...,yN)m^ 

Notice that coefficients Pj depend on time through $ 
even though the magnetic field is stationary. However, 
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FIG. 2. (Color Online). The coupling constants, Ji/Ji(0) 
as functions of the barrier height, a. The upper thick (blue) 
curve describes Ji , the lower (orange) curve corresponds to J 2 . 
With the upper most arrow, we show the limiting value of Ji 
for a —^ 00 . The insets show the system for the corresponding 
values of a. 


since this time dependence can be different from Ji{t) the 
probability of each ordering in the spin chain, i.e. |ai(t)p 
changes with time. An application of magnetic field can 
drive a transition to a spin segregated state where the 
spin down particles are not mixed with the spin up par¬ 
ticles. As was discussed in Ref. [s^, this is possible due 
to a high degeneracy of the spectrum such that even a 
tiny magnetic field gradient can drive such a transition. 
This can be utilized using the magnetic field in the form 
hix) = h^x, where &o is some constant parameter (there¬ 
fore Pj = A(t)/3j°^). By taking this constant to be large, 
such that /3j (0)/Jj(0) 1, we can have the initial state 

to be almost fully spin segregated or ’ferromagnetic’. By 
increasing the frequency of the external confinement we 
can drive the system from dominantly ’ferromagnetic’ to 
’antiferromagnetic’ states, since the Heisenberg Hamilto¬ 
nian in Eq. © is ’antiferromagnetic’. 

To conclude the presentation of the formalism, we con¬ 
sider a trapping potential where the quantum dynamics 
of a spin chain is altered without applying an external 
magnetic field. For this we use a potential schemati¬ 
cally shown in Fig. [1] where a shallow area with a time- 
dependent barrier in the center is surrounded by impen¬ 
etrable wells. We model this trap by an infinite square 
well potential, i.e. V{x, t) = af{t)6{x) for a; € [—1,1] and 
otherwise V{x, t) 00 . To give a spin chain time to react 
on the change of potential, we assume that f{t) varies sig¬ 
nificantly only on a time scale given hy gf. This assump¬ 
tion means that <h(t) changes almost adiabatically, which 
however does not imply adiabatic change of If (t) due to 
the degeneracy of the spectrum. Having this in mind, let 
us first assume that f{t) = 1 and study Ji for different a. 
Note that conservation of parity leads to Ji[a] = J 3 [q:]. 
So it is enough to study only the combinations Ji [a]/ Ji [0] 
and J 2 [a]/Ji[ 0 ] which are ^/-independent and are shown 
in Fig. [ 2 ] For a = 0, we have a pure infinite square 
well potential which requires Ji = J 2 . Positive values of 
a naturally descrease J 2 and increase Ji, such that for 
a —> 00 , we have J 2 IJiY^] —t 0 and Ji/Ji[0] —>■ 4/3. The 
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t/r 1/r 


FIG. 3. (Color Online.) Dynamics of probabilities for differ¬ 
ent configurations as a function of time for a four body system 
in an infinite square well potential with a time-independent 
barrier in the middle. Panel a) corresponds to A = 0.2, b) 
A — 0.5, c) A = 0.8 and d) A = 4. The solid (blue) 
shaded curve shows the change of J 2 /Ji[ 0 ] as a function of 
time. Dashed (orange) curve depicts ttf-f- configuration, dot¬ 
ted (red) - 4-4tt' These and all other configurations are shown 
in panels b), c) and d). Note that t4-4t ^-nd 4-tt4 configura¬ 
tions are equally probable. 

increase of Ji[q;]/Ji[ 0] is related to the increase of the 
density in one well by increasing the barrier. Note that 
this effect should be less visible for more particles. 

To illustrate the effect of this change of Ji[a]/Ji[0], we 
assume that for a —oo we prepare the system in the ttZ-'l 
configuration, see Fig. [TJ Now we open and close the 
barrier and investigate the evolution of the system during 
this cycle. We assume the following form of the variation: 
f{t) = (1 — \/sin(At/r)), where r = 7r/Ji[0] defines the 
natural time scale in the box in the absence of the barrier. 
Note that one cycle happens within the period 7r/(Ar). 
To supress the dynamics between different wells, we put 
a = 50 (see Fig. [5]). We present our hndings in Fig. O 
showing the probabilities of different configurations, i.e. 
joip for different A. Note that these probabilities after 
one cycle depend strongly on A which provides a way for 
state preparation. For example let us take a look at the 
case with A = 0.8 (panel c)). We see that for such a 
driving mode one ends up in 0.996% of all cases in 
which can be seen as the exchange of the pairs. Note, that 
if we had plotted the total density before and after the 
cycle for every configuration from Fig. |3l then we would 
have obtained the same result due to the adiabaticity 
of particle motion. Nevertheless, the spin configurations 
are profoundly different which highlights the separation 
of the spin and particle dynamics that we have derived. 
It is worthwhile noting that for even slower change of 
the potential with the time scale much larger than gf 
the dynamics in the both spin and charge sectors will be 
adiabatic. 
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Conclusions. - In this paper, we discuss a time depen¬ 
dent spin chain which is realized with strongly interact¬ 
ing atoms in a time-dependent confinement in ID. First, 
we outline a mapping onto a spin chain for one impurity 
in a Fermi sea of majority particles. This mapping can 
be trivially extended to more impurity particles or other 
multicomponent strongly coupled systems. Next, we use 
a time-dependent harmonic oscillator potential with a 
weak stationary magnetic field and an infinite square well 
potential with a time-dependent barrier to illustrate some 
basic properties of the spin dynamics in such systems. In 
particular, we show that in the former case by changing 
the trapping potential one can drive a system to a spin 
segregated state. For the latter case, we demonstrate the 
possibility of a state preparation and manipulation by 
proper changing the shape of the trapping potential. 


A major goal of cold atomic gas research is to reach the 
regime where quantum magnetism can be studied and a 
number of pmneering experiments have already been re¬ 
ported [SlrlSSj. In particular, the superexchange of two 
spins has been observed in Ref. 33| and it was shown that 
the lattice spin model limit of the Bose-Hubbard model 
could accurately describe the data. While lim¬ 


ited to strong cooupling and ID, the approach described 
here goes beyond those models as it can fully incorporate 
the shape of any (time-dependent) potential, circumvent¬ 
ing any need for making a lattice approximation. Our 
approach therefore ideally complements the lattice ap¬ 
proximation as a tool to simulate and study spin dynam¬ 
ics. This allows us to address the dynamical evolution of 
general A^-body exchanges in arbitrary potentials in the 
strongly interacting limit for both fermionic and bosonic 
atoms. Our theory may therefore be relevant for using 
exchange interactions to generate multiparticle entangle¬ 
ment 36|, l37| and building robust quantum gates 381 l39l| 
for use in quantum communication (4fll |. computation, 
and information (4lj |. 
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Supplemental Material for ’’Simulation of Time-dependent Heisenberg Models in ID” 


Here we outline the derivation of Eq. (5) of the main text. For this we first insert the wave function = (j)-\- -^f in 
the Schrodinger equation. Next we make a projection onto a specific ordering by integrating the equation with 
This procedure yields 

+ y^ixjmif) ■ ( 13 ) 

To proceed, we notice that Hxj^ = everywhere except at the points where the particles meet. We also notice 

that due to a non-smooth behaviour close to these points Hxj^ yields a Dirac delta function. These observations 
allow us to conclude that 

{x,^H\c|,)=^a,it){x,^^). ( 14 ) 

Next we turn our attention to the {xj^\H\f) term in Eq. (1131) which, for convenience, we rewrite as an integral over 
the configuration yi < y 2 < ■■■ < Uj-i < x < yj < ... < yN, 

/ OO ^OO POO poo poo 

dyi / dy 2 ... / dx dyj... / dyNdx...dyN^*Hf. (15) 

-OO Jyi Jyi-i JyN-i 


Our next steps are two integrations by parts. This will yield some boundary terms and the integral with / and $ 
exchanged. Notice that there will be two types of boundary terms: i) with yi = yk, and ii) with x = yj-i or x = yj. 
The former terms vanish due to the fermionic nature of the majority particles. The latter, however, should be properly 
taken into account. 


ixjmif) 

m 


-lim 

2 £->-0 


yi<y2<...<yj-i<x<yj<...<yN 


+6ix-yj-i -e) 


da;d 2 /i...dyAr( [(5(a; - yj -be) - 6{x - y^-i - e)] ( 


df 


-f-. 


d^* 


dyj-i dyj-i 




( 16 ) 















6 


where all but the last terms under the integral sign are the boundary terms. Thus, we have 

{Xj‘^\H\f) = -i{fXj\-^)* + boundary terms . 

Next, we collect the expressions just derived and write the equation for aj{t) 

. duj {t) 1 boundary terms 


dt 


9f 


( 17 ) 


(18) 


where we used that = 0 by construction. It is important to notice that Eq. (fTSl) is general and does not rely 

on the assumption that 1/gf is small. However, as we show below this assumption makes the derived expression very 
useful. Let us now focus on the boundary terms. 


TV! 

boundary terms =- — lirn 


dxd2/i...d?/7v/x 


5[x-yj+t) 


yi<V2<...<yj-l<x<yj<...<yN 

d^* d^* 


dyj 


dx 


+ 5[x-yj-i -e) 


d^* 


dx dyj-i 


N\ 


dxd?/i...d?/Ar/( (5(a; - yj) - 5[x - yj-i)\-^. 




(19) 


J yi<y2<...<yj-i<x<yj<...<yN 

To obtain / at the points where the particles meet, we use the boundary conditions from the main text. This yields 


= yi) 






1 

‘^9f 



9yi) x=yf 


V dx dyi j 


( 20 ) 


Now we make the assumption that l/ff/ 0, which implies that f{x = yi) ~ (a^+i — ■ Inserting this result 

in Eqs. m and (HU), we arrive at the desired expression. It should be noted that in our derivations we assume that 
gf sets the largest energy scale of the problem. Therefore, if the change of the trap is such that gf becomes of the 
order of the other energy scales then the treatment above is not valid, namely we cannot neglect the second term on 
the right-hand-side of Eq. (l20)) . Eor example this can happen if we increase the density of the system, by squeezing 
the trap, which necessarily increases the kinetic energy. Another instance is a periodic driving in the parametric 
resonance region, which pumps in energy in the system. 
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